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We discuss an efficient scheme for obtaining spin-polarized quasi-particle excitation energies within 
the general framework of the density functional theory (DFT). Our approach is to correct the DFT 
eigenvalues via the electrostatic energy of a majority or minority spin electron resulting from its 
interaction with the associated exchange and correlation holes by using appropriate spin-resolved 
pair correlation functions. A version of the method for treating systems with localized orbitals, 
including the case of partially filled metallic bands, is considered. Illustrative results on Cu are 
presented. 
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As spectroscopic data using a variety of experimental 
techniques is becoming available at an ever increasing 
pace, there is great need to develop efficient methods for 
calculating electronic excitation energies in wide classes 
of materials. The GW approximation (GWA) has been 
the traditional route for addressing this problem 1 . The 
GWA however is computationally quite intensive, which 
makes the routine application of the GWA in electronic 
structure computations difficult and for this reason GWA 
results have to date been limited largely to relatively sim- 
ple systems. 



With this motivation, in a previous stud}^, we have 
attempted to develop a simplified scheme for obtaining 
quasi-particle energies by correcting the eigenvalues e n 
given by the density functional theory (DFT)— with a 
state-dependent correction A„ associated with the self- 
energy Yj xc of the exchange-correlation hole surround- 
ing each electron. The results are similar to those ob- 
tained within other extensions of the DFT 4 i 5 i 6 i 7 i 8 iV°, 
which have been shown to describe the energy gaps in 
semiconductors with an accuracy comparable to that of 
the GWA. The purpose of the present work is to gener- 
alize the unpolarized case considered in Ref. to include 
spin-polarization effects. We also address the treatment 
of self-interaction corrections in systems with localized d 
and / orbitals and suggest a scheme for handling metallic 
systems with partially filled bands of d or / character. 



This article is organized as follows. The introduction 
is followed in Section II by a brief explanation of the rel- 
evant pair-correlation functions and their role in evaluat- 
ing self-energies. The spin-resolution of the exchange and 
correlation energies is developed in Section III. A some- 
what more general formulation of our scheme is taken up 
in Section IV which also presents illustrative results in 
Cu. Finally, Section V makes a few concluding remarks. 



We consider the elements 

7r(n, r 2 ) = N(N - 1) J |*(n, r N )\ 2 dr 3 ...dr N , (1) 

of the two-particle density matrix, where N is the total 
number of electrons and VP is the ground state wavefunc- 
tion of the many body electronic system. The pair cor- 
relation function g(ri, r 2 ) is then defined via the relation 



7r(ri,r 2 ) = n(r 1 )n(r 2 )g(r 1 ,r 2 ) 

= n(n)n(r 2 )(l + C(r 1 ,r 2 )), 



(2) 



where n(r) denotes the electron density and the corre- 
lation factor C(ri,r 2 ) is defined by the second equality. 
For an uncorrelated system <?(ri, r 2 ) — > 1, or equivalently, 
C(ri,r 2 )-»0. 

Two effects must now be considered, namely, exchange 
and correlation. The antisymmetry of the many electron 
wavefunction ^ prevents electrons of the same spin to 
come too close to one another and creates an 'exchange 
hole' around each electron. The associated exchange en- 
ergy is: 



where 



E r 



e*(r) 



dr n(r)e x (r) 



I f dr , n(r')C x (r,r') 



(3) 



(4) 



is the exchange energy per particle and C x is the ex- 
change contribution to the pair correlation C(r, r') of 
Eq. [H The factor of 1/2 accounts for double counting 
of the electron-electron interaction. Thus, the exchange 
energy per particle e x (r) can readily be interpreted as the 
Coulomb interaction of an electron at r with its exchange 
hole. For the homogeneous electron gas (HEG) of den- 
sity n, C x and e x have well-known analytical expressions 
C x EG (n) and e^ EG (n) and, in particular, the familiar 
Hartree-Fock-Slater (HFS) potential is given bjiii 



v s {v)=2e^{n{v)) 



(5) 
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Here the factor of 2 arises from the fact that variations 
with respect to a trial wavefunction of the form of a Slater 
determinant yield a potential term in one-particle equa- 
tions which is twice as large as the corresponding en- 
ergy term e x . The DFT within the local density approx- 
imation (LDA) provides a remarkably similar scheme in 
which the electron density n(r) is treated as the varia- 
tional parameter—. The exchange hole is still described 
by the pair correlation function between like spin elec- 
trons, but the exchange potential 



v x (r) 



SE X 
Sn(r) 



(6) 



differs from the HFS potential by a factor of 2/3. 

Coulomb repulsion between electrons is an additional 
effect responsible for electron motions to become corre- 
lated. As a result, there also appears a 'correlation hole' 
around each electron, which is dominated by electrons 
of opposite spin, since the like spin electrons are already 
excluded via the exchange hole. The correlation hole 
is generally less deep compared to the exchange hole, 
especially at high electron densities. Note that the ex- 
change hole goes to zero at the origin and excludes a to- 
tal of one electron worth of charge in order to correct for 
self-interaction terms. In contrast, the correlation hole 
integrates to zero and involves only a redistribution of 
charge^. 

A number of intuitive schemes for treating the ex- 
change and correlation holes have been proposed, where 
the pair correlation functions are determined from a 
Schrodinger equation for an electron pair with appro- 
priate boundary conditionsii2ii£ii£iiLi2i Such models are 
useful for generalizing the LDA by improving the descrip- 
tion of exchange and correlation in inhomogeneous sys- 
tems. Note, however, that in going beyond the LDA, 
a distinction should be made between theories which at- 
tempt to find better energy functionals but still lie within 
the framework of the DFT, and those theories which fo- 
cus on the self-energy £ xc in order to obtain better quasi- 
particle energies. The present scheme belongs to the lat- 
ter category where the excitation energies are given by 



E,, 



(7) 



Here e n is the eigenvalue corresponding to the Kohn- 
Sham orbital ip n , and A„ is a state-dependent correction 
associated with the self-energy E xc . By using a particular 
ansatz for Y, xc , we have previously proposed 2 



A„ = J d 3 r (2e xc (r) 



v xc {r)) \Mr)f 



where 



v xc (r) 
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Sn(r) 

is the Kohn-Sham potential and 



E xc = / dr n(r)e xc (r) 



(8) 



(9) 



(10) 



is the exchange-correlation energy. The exchange- 
correlation energy per particle e xc can be expressed in 
terms of the pair correlation functions as 



=(r) = 



1 



dr' 



v^y), (id 



where A is the coupling constant from the Hellmann- 
Feynman theorem. Within the LDA, e xc (r) is a well- 
known function of the local density n(r) obtained from 
homogeneous electron gas results^. 



III. SPIN RESOLVED QUASI-PARTICLE 
CORRECTION 



In the spin resolved case, the prescription for the cor- 
rection to the excitation energy given by Eq. |S| becomes 



A£ = J d 3 r (2e xc (r) 



where 



< c (r) 



< c (r)) |<(r)| 2 



8 E xc 



8n a (r) 



(12) 



(13) 



is the spin-dependent DFT potential and n a the spin- 
resolved electron density. As in Eq. 11 above, the 
exchange-correlation energy per particle e xc (r) is given 
by the Coulomb interaction of an electron of spin a at r 
with its exchange-correlation hole: 



2 



dr' 



dX- 



,(r')Cr ( r > r ') 



(14) 



in terms of the spin-resolved correlation functions 
(r, r'). It is convenient to split e xc (r) as 



4 c (r)=4(r) + <(r), 



(15) 



where e x and e° give the exchange and correlation com- 
ponents of e xc (r), respectively. The spin-resolution of 
e x is straightforward because the exchange hole only in- 
volves like spin electrons. On the other hand, the spin- 
resolution of e° is more subtle since here both like and 
unlike spin electrons contribute. The energy per particle 
for the spin polarized HEG depends only on the density 
parameter r s and the local spin polarization 



C = (n T - ni)/n 



(16) 



In this connection, the spin-resolution of the cor- 
relation energy presented recently by Gori-Giorgi and 
Perdewi^for the HEG may usefully be employed. Ref.ll9l 
defines spin resolved correlation energies 



C (r s ,C) = 27r-^ 



T ,C°> a (r a ,C,u) udu , (17) 
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FIG. 1: Characteristic behavior of the spin-resolved 
exchange-correlation energies ej c and e^ c and the correspond- 
ing spin only contributions ej and q as a function of the spin 
polarization £ for r s = 2 in a HEG based on Eqs. 15-18. See 
legend for the meaning of various line types. 



FIG. 2: Spin- resolved correlation energies e^(r s ,Q (see Eq. 
18) in a HEG are shown in surface plots as a function of the 
spin polarization £ and the density parameter r a . Upper and 
lower surfaces refer to the majority ( j) and minority (J,) spins, 
respectively. 



by using correlation holes n a iC° a (r s , (, |r — r'|), where 
the pair correlation function C° ,a has been averaged over 
the coupling strength A. Here we have made the depen- 
dence of various quantities on r s and C explicit. Thus 
the interaction energy e° of an electron of spin er with its 
correlation hole can be resolved as 



n 



(18) 



Parameterized expressions for e° ,a and e£' a based on 
the HEG are given in Ref. [H 

Figs, n and present illustrative results concerning 
the r s and £ dependencies of the exchange and correlation 
energies in the HEG. Fig.[T]shows the behavior of ej c as a 
function of £ for r s — 2. The correlation contribution t a c 
here has been obtained from Eq. El The spin-splitting 
of e xc is seen from Fig. Q] to increase steadily with £. 
The exchange contribution at a given density depends 
linearly on (1 ± C) 1 ^ 3 so that when £ is small t% displays 
a linear behavior with £. Notably, t a x for minority spin 
electrons is seen in Fig. Q] (thin dashed line) to vanish in 
the ferromagnetic limit (£ = 1) reflecting the fact that a 
minority spin electron has no other electrons to exchange 
with in this case. In contrast, the correlation energy per 
particle of the minority spin electrons does not vanish 
in the limit £ = 1 due to the contribution of opposite 
spin electrons in Eq. El The effect of correlations is to 
reduce the £ dependency of e£. As a result, in the linear 
region, the slopes of ej c curves are seen to be smaller and 
magnetic splittings are significantly reduced. 

In order to further highlight the role of correlations, 
we show ej (upper surface) and t\ (lower surface) as 
a function of the parameters £ and r s in Fig. [3 It is 



clear immediately that correlations yield a less attractive 
hole for the majority spins compared to that for the mi- 
nority spins. The reason is that correlation effects are 
dominated by the contribution of opposite spin electrons 
which are not prevented by the Pauli principle to come 
close to each other. The splitting of the two surfaces in- 
creases with increasing spin polarization £ but decreases 
with increasing r s . At C = the two surfaces collapse on 
to the well known paramagnetic curve, which scales asli 



0.127 
e c = — Hy 



(19) 



for typical metallic densities (r s = 1 — 6). 



IV. A MORE GENERAL FORMULATION 

This section presents a framework for generalizing 
Eq. [7| This formulation also allows us to connect our 
approach to some related methods in the literature. In 
the interest of notational simplicity, we will omit from 
hereon the obvious spin dependencies and the associated 
spin-indices. 

We recall that the standard definition of the the self- 
energy operator S xc is: 



Gn 



G~ 



(20) 



where G is the exact single-particle Green's function and 
Go is the reference Green's function obtained within the 
Hartree approximation. Thus T, xc in principle contains 
all quasiparticle effects beyond the Hartree approxima- 
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tion. We now write T, xc (suppressing spin-indices) as 



^xc 



V X c + A„|t/i„ >< tp n \ 



(21) 



where v xc is the spin-dependent exchange-correlation po- 
tential associated with the spin-dependent Kohn-Sham 
orbitals 4>n~ The state-dependent correction A n is more 
general than the expression of Eq. ^] and it can be the 
basis for curing the unphysical interaction of an electron 
with itself which occurs in the DFT. In fact, the right side 
of Eg. 1 121 attempts to correct for such self- interactions by 
using the Coulomb potential of the exchange-correlation 
hole, but does so only approximately since it assumes 
every electron to possess the same pair-correlation func- 
tion. In systems with localized orbitals, it makes sense to 
subtract out the self-interaction for each occupied orbital 
explicitly by usingSi 



A„ = 



J d 3 r (- 



■MO-O)) \Mr)\ 2 



(22) 



where u n (r) and u" c (r) are the electrostatic and the 
exchange-correlation potential, respectively, created by 
the charge density |V> n (i")| 2 - The self- interaction cor- 
rection (SIC) given by Eq. has been proposed by 
Perdew and Zunger— for atomic systems and possesses 
a straightforward spin-dependence. For atoms, SIC low- 
ers the LDA eignevalues and improves agreement with 
measurements^. Note, however, that in extended sys- 
tems SIC encounters conceptual problems due to the de- 
localized character of the Bloch functions^, In fact, the 
SIC correction turns to be zero for plane-waves. A route 
around this problem is to renormalize the wavefunction 
within the unit cell and perform the correction within 
the atomic spheres i24 

The so-called LDA+U can be considered a simplified 
version of the SIC where the correction A^ acts only on 
localized d-oi /-states, often splitting the metallic LDA 
bands into upper and lower Hubbard bands descriptive 
of Mott-Hubbard insulators25i2Si2i. However, problems 
with both SIC and the LDA+U arise for systems with 
partially filled 3d shells which are metallic, as for exam- 
ple in the transition metalsi2£ In order to handle cases 
where a conduction band crosses the Fermi energy Ep, 
we observe first that here the DFT potential v xc generally 
provides a reasonable approximation to the self-energy at 
the Fermi surface (FS)2& and we may impose the condi- 
tion 



^xc\FS — v xc 

on the self-energy. We now write A n as 
A„ = A + A„ , 



(23) 



(24) 



where A is defined to be the average correction on the 
FS and A„ is a state-dependent contribution. It is clear 
that the condition of Eq. [23 can be satisfied by setting 
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FIG. 3: State-dependent correction A n based on Eqs. 1241 
and 1221 in Cu as a function of energy with zero denoting the 
Fermi energy. Each dot gives the value of A n for a particular 
state. Total number of states shown is 6 x 505 corresponding 
to 6 valence bands at each of the 505 k-points computed in 
the l/48th irreducible portion of the Brillouin zone of Cu as 
discussed in the text. 



the average shift A to zero and requiring that A n goes 
to zero at the Ep, which leads us to suggest the form 



K = A n 



e n — Ep 
Eb — Ep 



(25) 



in metallic systems. Here Eb denotes the bottom of 
the conduction band. The linear interpolation given by 
Eq. 1251 is equivalent to a band renormalization near the 
Fermi level given by 



U = (1 + A n )£„ 



where 



A ra — 



Eb — Ep 



(26) 



(27) 



and £„ = (e n — Ep). Such a renormalization due to self- 
energy effects is often invoked phenomenologically within 
the framework of Fermi liquid theories^ 

Fig. |31 shows illustrative results in Cu where the elec- 
tronic structure involves states of a mixed s — p and d 
character. For this purpose, the standard LDA-based en- 
ergy bands and Bloch wavefunctions were computed for 
the 6 valence bands of fee Cu on a uniform 505 k-point 
mesh in the irreducible l/48th of the Brillouin zone. The 
correction A„ of Eas. 1221 and 1241 was then computed for 
each of the 505x6 states^. The resulting values of A n 
after the average A is subtracted are shown in the 'dot- 
plot' of Fig. |21 A„ is seen to be positive and fairly con- 
stant at low band energies (—0.7 to —0.5 Ry), where the 
states mainly possess s character. At intermediate ener- 
gies (—0.4 to —0.2 Ry), where the d admixture increases, 
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A„ decreases rapidly as self-interaction effects become 
stronger. Closer to the Fermi energy (—0.1 to 0.0 Ry), 
the s — p character increases once again and this is re- 
flected in the upturn in A„ values seen in Fig. |31 Using 
Ep — Eb = 0.691 Ry, we find that near the FS, typical 
A n values are « 10 %, which is consistent with the value 
of the renormalization parameter A f=a 8 % adduced from 
fitting experimental photoemission results2ii22i2i. 

In ferromagnetic metals such as Fe and Ni the present 
scheme leaves the band structure at the Fermi energy 
and hence the majority and minority spin FS's and the 
corresponding spin magnetic moments unchanged due to 
the constraint of Eq. 1231 However, energy bands away 
from the Fermi level arc modified. In particular, the top- 
most occupied d states in 3d and 4c£ transition metals, 
which are more localized within the atomic sphere than 
the d states toward the bottom of the band, will expe- 
rience larger corrections (lowering). As a result, d-band 
widths would become narrower giving trends consistent 
with photoemission experiments^. 



for the homogeneous electron gas. The approach is sim- 
ilar to that followed in our earlier treatment of the un- 
polarized case^, and should yield quasi-particle energies 
with an accuracy comparable to that of the more de- 
manding GWA. Notably, our scheme does not involve 
the limitations of the standard implementations of the 
GWA in which the screened interaction W is assumed 
to be spin-independent. We also show how our scheme 
can be viewed in a somewhat more general context, so 
that quasi-particle energies can be corrected on a state- 
by-state basis for self-interaction terms in the DFT along 
the lines of the familiar SIC scheme^. The question of 
applying self-interaction correction to a metallic band is 
considered and an interpolation formula to handle such 
a case is suggested. Finally, illustrative results in Cu are 
presented and the computed renormalization of the DFT 
eignvalues of the conduction band in Cu is found to be in 
reasonable accord with the corresponding experimental 
results. 



SUMMARY AND CONCLUSIONS 
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We discuss how DFT eignevalues can be corrected via a 
spin-dependent self-energy Yf xc to obtain improved quasi- 
particle energies. For this purpose, Yf xc is evaluated as 
the Coulomb energy associated with the exchange and 
correlation holes surrounding each electron in terms of 
the appropriate spin-resolved pair correlation functions 
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